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/\ ' Abstract 



We investigate the possibility that both the baryon asymmetry of the universe and 
the observed cold dark matter density are generated by decays of a heavy scalar field 
which dominates the universe before nucleosynthesis. Since baryons and cold dark matter 
have common origin, this mechanism yields a natural explanation of the similarity of 
the corresponding energy densities. The cosmological moduli and gravitino problems are 
avoided. 



1 Introduction 

Supersymmetry is a well-motivated framework for TeV scale physics and it has several 
conceptually nice features. Apart from being the unique extension of the space-time 
symmetry it provides us with the most compelling scenarios of gauge coupling unification. 
Not only for short-distance physics, supersymmetry provides interesting possibilities for 
the explanation of the structure of the Universe. The supersymmetric Standard Model 
contains a natural candidate for cold dark matter of the Universe, and also the fact 
that quadratic divergences are absent allows us to naturally postpone the cut-off of the 
theory as high as the Planck scale. It makes us possible to calculate high-scale or high- 
temperature phenomena such as inflation and reheating in a reliable framework. 

However, a closer inspection of cosmological scenarios for dark matter production 
and baryogenesis in supersymmetric models reveals that there is often an inconsistency 
in the underlying assumptions. For example, the most popular scenarios for dark matter 
and for baryogenesis are known to be incompatible. It is widely accepted that the relic 
density of the neutralino from thermal decoupling naturally explains the amount of cold 
dark matter, and there is a good explanation of the baryon component of the Universe by 
the thermal leptogenesis scenario ]l|]. The thermal neutralino dark matter realizes in the 
gravity mediation scenario, i.e., the gravitino mass m^/2 ~ 100 GeV, whereas the thermal 
leptogenesis needs a high reheating temperature after inflation, Tr > 10 9 GeV [Q, ||, |4|. 
Considering the thermal production of the gravitinos and their decays, the bound from 
the big bang nucleosynthesis (BBN) excludes such a possibility (see e.g. || ||, |7[ ||, |9|). 
It is possible to avoid the constraints from the BBN by assuming a very heavy (rn.3/2 •> 
100 TeV) or a stable gravitino, but the most popular explanation of dark matter via 
the thermal relic neutralinos will be lost in such cases since the lightest supersymmetric 
particle will be the lUino (l^, |ll[ (or Higgsino (T^| ) , which typically leads to a too small 
abundance, or the gravitino. 

Not only that, for m-3/2 ~ 100 GeV, there is a serious cosmological moduli (or Polonyi) 
problem [13]. It is possible to give a large mass to the Polonyi field by assuming a non- 



trivial interactions in order to avoid the late time decay [14], but in such a case the decay 



into gravitinos causes a disaster [15, 16]. There is also a rather disturbing bound on the 
reheating temperature coming from modular cosmology [Q|. 

It is clear from the above discussion that we need to abandon some of the "standard" 
assumptions. We propose in this paper an extreme but cosmologically safe scenario in 
which radiation, dark matter and baryon asymmetry are all generated right before the 
BBN era so that there is no constraint from the high-temperature gravitino production, 



moduli dynamics etc. (See [18, 19, £(], 21, |22[ for earlier works). We find that there is a 
consistent framework to realize this scenario by a late-time decay of a scalar condensate 
(j) in supersymmetric models. We assume that coherent oscillations of a scalar field first 
dominate the energy density of the universe. The <j) field then decays later producing a 
large amount of radiation and all the dangerous relics (e.g. gravitinos) get diluted to a 



negligible density. The standard BBN takes place after <$ has decayed at sufficiently high 
temperature. The dark matter is produced by the <p decay and the correct abundance is 
obtained in the H^ino/Higgsino LSP scenario due to pair annihilation p3| ] which is less 
effective than in the so-called thermal scenarios. The baryon asymmetry is also generated 
in (j) decays. The decay channel of a <f> field into two gravitinos is forbidden by -R-parity, 
and we assume that the decay into a gravitino and the fermionic superpartner of eft 
((ft) is kinematically forbidden. Interestingly, by unifying the origin of dark matter and 
baryon asymmetry, this scenario explains the one of the puzzling issues of our Universe 
that the energy densities of dark matter and baryon asymmetry are close to each other, 

The paper is organized as follows: in Section || we present the baryogenesis mechanism, 
and show that it is indeed possible to generate the observed amount of baryon asymmetry 
while satisfying the BBN constraints on the reheating temperature after the (f> decay. In 
Section 0, we discuss the abundance of the dark matter from the <f> decay. In particular, 
in our preferred scenario the ratio f2cDM/^& ~ 5 implies a large gravitino mass, m^/2 ~ 
100 TeV which fits nicely to the lyino/Higgsino LSP scenario. 

2 Baryogenesis 
2.1 Basic idea 

We consider a chiral superfield <£ = ((f), (f>, F$) which couples to the matter fields via a 
higher-dimensional term in the superpotential f20[|, 



W D -&UDD, (1) 

M K ' 

with U = (u c ,u c ,F u ) and D = (d c , d c , F^) denoting the up- and down-type quark su- 
perfields. Here, we suppressed color and generation indices, and absorbed dimension- 
less couplings into M which will be taken to be of the of order of the Planck scale, 
M ~ Mp = 2.44 x 10 18 GeV, unless stated otherwise. Due to this operator <fi effectively 
carries baryon number (+1). 

Let us now define the <j> number asymmetry 



q§ is given by the difference between the number densities n^ and n^* of particles (f) and 
antiparticles (f>*. qs can be interpreted as angular momentum of the <j> field rotating in 



the complex plane [29]. 



The scenario we shall describe in the following consist of the following sequence of 
steps: first, a positive q<h is generated. Then there is an era of coherent </> oscillations where 
a significant fraction of the energy density of the universe is carried by these oscillations, 



1 See, for example, [E4l Eq, E6l E71 Ell E8[ for earlier attempts to explain the similarity: f^ ~ Q 



CDM- 



and q^ is conserved. Finally, <f> number is converted into the baryon number by its decay, 
and it reheats the universe up to 0(100 MeV) consistently with nucleosynthesis. 

Before we describe the mechanism in detail, let us briefly explain the main differences 



to Ref. [20]. For the mechanism to work, one has make sure that dangerous ^-number 
violating interaction terms are absent or sufficiently suppressed. One can forbid these 
terms by imposing a symmetry. In Ref. |2(|, a model with a Z4R symmetry is presented, 
which ensures the ^-number conservation at a sufficient level and also prevents (j) from 
dominating the universe. Below, we will consider a different model with an anomaly-free 
Z9 discrete baryon symmetry in addition to the usual iJ-parity (see Table ||), and will, 
as already stated, assume that that eft dominates the universe at an early epoch. 

In order to obtain the baryon asymmetry before BBN, in Ref. p3] enhanced couplings 
of <p to the baryons are assumed such that the eft lifetime is short enough. Moreover, in 
the case of the <f> mass of the order 100 GeV, the universe would be always matter (</> 
or dark matter) dominated once (j) dominates the universe at an early time. This is 
not compatible with the requirement of successful BBN and, therefore, a model without 
(early) <f> domination is constructed there. 

The situation is, however, different if <f> is heavy. As the temperature after the <f> decay 
is higher than the BBN temperature, cj> is allowed to dominate the energy density of the 
universe. This high temperature after <f> decay also plays a crucial role in the generation 
of cold dark matter (cf . Section ||) . 

2.2 <fi evolution 

Let us start by considering the dynamics of the (f) field. The evolution of (j) is described 
by its equation of motion, 

dV 

+ ( 3 # + W+ = 0, (3) 

where V = V(<j>, 0*, . . . ) denotes the scalar potential, H the Hubble rate and 1^ the <j> 
decay rate. Eq. (|3|) translates into an equation of motion for q^, 

( dV dV\ 

fc +Sff »- -'(*5r*v)- (4) 

Hence, a non-vanishing right-hand side of (||) can be used for the first step, i.e. to create 
non-zero q^ dynamically. Before explaining this in detail, recall that we need also to 
satisfy the condition of (p number conservation in the stage of <fi oscillation. This means 
that in the eft oscillation era the </> number violating terms have to be absent (or sufficiently 
suppressed). The most dangerous term of this type is fi 2 (p 2 + h.c. 

In order to enforce the absence of those dangerous terms, we impose a discrete Zg 
symmetry which is an anomaly free subgroup of baryon number symmetry [BH, 31, 32, 
|33| 34 1 . The charge assignment is listed in Tab. |l|, where <I> is introduced in order to give 



a mass term for the fermionic superpartner of <f> without introducing fj?(ft 2 term in the 



field 


Q 


u 


D L E N H H $ 


$ 


.R-parity 


+1 


-1 


-10 0+3 

---- + +- 


-3 



Table 1: Charge assignments under Zg and i?-parity. We denote the MSSM superfields 
according to Q = (q, q, F q ) with q representing the quark doublets etc. 

Lagrangian.fj With this choice (j) = can always be a minimum of the potential which is 
necessary to preserve the imparity. 

The symmetry does, however, allow for a <3? 6 term in the Kahler potential and in 
the superpotential. In the following we discuss the two cases in which the $ 6 term in 
superpotential is absent and present. 

Case A: no $ 6 term in the superpotential 

In the case where there is a $ 6 term only in the Kahler potential, the potential for the 
cf) field is given by 

ml /2 M 2 F{\^\ 2 /M 2 -- 

+ higher-order terms , (5) 



V = ml\4>\ 2 + ml /2 M 2 F(\4>\ 2 /M 2 



« -MZ- 4> 6 + h.c. 



M 4 

where k is expected to be order one. The gravitino mass parameter m 3 / 2 represents the 
supersymmetry breaking scale. The function F(x) is a general (polynomial) function. 

The presence of the </> 6 term in (||) can lead to a dynamical generation of q^ ^ as 
follows: for H S> ms, 4> oscillation is negligible and we can treat (f> to be constant. We 
can integrate the equation 

2 
777 

qt + ZHqt = -^Im[ K 6 ] , (6) 

so that 

2 
777o /o n 

qAt = 777T ) ~ \k\ 7 d>f„: , (7) 

W 4, ) 11 2m0 M 4 ^ m ' v ; 

where 0i n ; is the initial amplitude of 4> after inflation which is generically O(M) with the 
potential in Eq. (||) (with rni, 2 replaced by 0(H 2 )). The number density of <fi and (f>* 
particles is given by p^/m^ where p^ ~ m\ \4>\ 2 + \4>\ 2 . Since the oscillation of (ft starts 
with amplitude (ftim, we obtain for the dimensionless (ft asymmetry 

«l f^V • (8) 



n</, + Tip \ 777^, 



2 The analysis will remain unchanged when we include the dynamics of the (ft field. Although a possible mass 
mixing term, m 2 (j)(f> + h.c, will distribute the baryon number to the </> field, Eq. (Eh will be the same once we 
include the (ft field in q^ in Eq. (0) . The decay of (ft (or more precisely the other mass eigenstate) happens about 
the same time as the (ft decay provided the mixing is order one. 



Here we have taken <p[ n { ~ M. We have checked that this expression yields roughly the 
correct order of magnitude for e (as long as 0i n j is comparable to M) by solving the 
equation of motion numerically. 

The resulting cf) asymmetry stays constant after <fi starts to oscillate because of the 
r.h.s. of Eq. (g) becomes numerically irrelevant when the amplitude drops (far) below 
M. As a consequence, n?q<f, (with R being the scale factor) is approximately conserved 
during the eft oscillation era, until eft decays. 

Case B: <fr 6 term in the superpotential 

The case with $ 6 term in the superpotential is qualitatively the same as case A. The 
potential of the <f> field in this case is given by 

r'2\ 



V = ml\cP\ 2 + ml /2 M 2 F(\cP\ 2 /M 2 



^3/2^ ,,.„!, „„H 10 



K M^ + hX 



+ k — + higher-order terms , (9) 

M 6 



where k' and k" are O(l) coefficients. In the early universe, m 3 / 2 which represents the 
SUSY breaking effect is replaced by the Hubble rate H. Now the minimum of the potential 
is generically 4> ~ {HM 3 ) 1 ' 4 which sets the initial amplitude of the 4> oscillation to be 
0ini ~ (rn^M s ) ' since the (j> oscillations start when H ~ m&. 

The equation for the evolution of the (p number asymmetry q^ in Eq. (||) is 

q^ + SHq^ = I^Im^'/] . (10) 

With the value of at the minimum of the potential, <fi ~ (HM 3 ) 1 ' 4 , before cf) oscillation, 
we find 

„ |^| m 3/2 0f„i nu 

m<f, M 6 
Therefore, the asymmetry factor e is estimated to be 

e~W\^, (12) 

which is larger than the case without the $ 6 term in the superpotential if 7713/2 ^ m^ 
and k' ~ 1. In conclusion, the presence of the <3? 6 term leads only to a quantitatively 
different result. 



2.3 Baryogenesis via <fi decay 

So far we have seen that, due to the presence of higher-order terms, a <p number asym- 
metry is induced which is conserved in the regime of <j> oscillations until eft decays. Let us 
now consider the conversion of <f> number to baryon number of the universe through the 
decay arising from the coupling ffih, <j> — > qqq. The corresponding decay rate is given by 

r> = ^, (13) 



where £ is obtained by a standard calculation. In the simplest case where all couplings 
of 4> to the quark superfields equal one, we obtain £ = 27/(2567r 3 ) ~ 3 x 10 -3 . 

The temperature Td of the thermal bath after eft decay is calculated by equating 
Hubble rate H and Fa, 

Td » ooMev ( « V' 2 ( a f 2 (Ml) . (14 ) 

V 10 " V Vl500TeV/ V M / 



Since has a large hadronic branching fraction, T^ has to fulfill [35| T^ > 4MeV. The 
corresponding lower bound on the eft mass is 

/ 10 -2\ 1/3 / M \ 2/3 
m ,>150Tev(— ) (-) , (15) 



where we have exploited that y^Sr — 1 for temperatures of the order MeV. 

Assuming that <p dominates the energy density of the universe before its decay, the 
number density of (f> just before decay is obtained by 

7T 2 

viia (rty + n r ) ~ — g* T d . (16) 

Using the relation between baryon asymmetry and number, n^ = q^ = e(n§ + n^. ), we 
can estimate the baryon asymmetry as 

3 3 (17) 



S 4 7710 



1/2 / rn 3/2 n2 / m x-3/2 /Mp 



10 10 ' N [w^ ) IsOTevJ 1 1500 TeV J (,lf ) ' (A) 

1U |K| VlO-V V2TevJ ll500TevJ ^ M ) ' (j 

Cases A and B correspond to the model without and with the <I> 6 term in the superpo- 
tential, respectively. By comparing the observed value (rib/s) = (8.7 ± 0.3) x 1CP 11 



1 36], the lower bound on ttu in Eq. ( |15[) requires either enhanced SUSY breaking terms 
|k| S> 1 or the large gravitino mass 7713/2 > lOTeV for case A whereas 1713/2 ~ ITeV is 
possible in case B (assuming M ~ Mp). 

As we shall see in the next section, the reference value of nriA ~ 1500 TeV is motivated 
by considerations on dark matter. Notice that the required size of the SUSY breaking 
scale 7713/2 from baryogenesis is indeed at a favorable value for low energy phenomenology. 
In particular, in case A, 7713/2 ~ 50 TeV indicates that loop corrections to soft masses 



(anomaly mediation) are important. In many models with such situation, the lUino [11 



[T0| ] or Higgsino |37], |38|, 39] becomes the LSP and that is indeed consistent with the 
discussion below. For case B, the Wino or Higgsino LSP is realized by either assuming 
small value or phase of k' such that gravitino mass is enhanced or simply assuming the 
Wino or Higgsino LSP in the scenarios of the gravity mediation type by relaxing the 
universality of the gaugino and/or scalar masses. 



3 Dark matter from (f> decay 

Since every (ft decay produces (at least) one superpartner, Ref. |2(| concludes that the 
number density of LSPs exceeds the one of baryons, tilsp > n^. However, tt-lsp is modi- 
fied by LSP pair annihilation processes in a heavy (ft scenario. These processes are effective 
as long as the corresponding rate exceeds the Hubble rate. In the MSSM, the dark matter 
candidates which pair annihilate strongly are the VKino and the Higgsino. Both particles 
have large annihilation cross sections through weak interaction, and thus the thermal 
abundance cannot explain the energy density of the dark matter. On the other hand, 
non-thermal production from (ft decays renders the H^ino/Higgsino a viable dark matter 
candidate. 

In order to find out to what extent the LSPs annihilate, one describes the evolution 
of number densities of (ft quanta and LSPs, n^ and n x , and energy density of the thermal 



bath, p ra d, by Boltzmann equations 22], 
dn, 



dt 



+ 3Hn (f) = -T^n^ , (18a) 



dn 



— -^ + 3Hn x = i/ LS p T^ n^ - (a v) n\ , (18b) 

— ^g-+4:Hp rad = (m^-^spm^T^n^ + m^av)^ . (18c) 

^lsp denotes the number of LSPs produced by a eft decay. The Boltzmann equations can 



be integrated, and the relic density of % can be approximated by |4C] 



^ ~ (A(av)M P T d )- 1 , (19) 

s 

as long as the LSPs are not equilibrated, i.e. Tj < m x /30. The relic x abundance is then 

o 7,2 „ Q1 p-5xlQ- 3 \ /1Q- 2 \ 1/2 / m x y / m<f) yt/2 f M_\ 
x ' V m 2 x (av) ) V £ ) VIOOGeV/ V1500TeV/ V m p/ 

(20) 

The thermal average of the annihilation cross section is typically (av) ~ 10~ 3 /m 2 for the 
particles which have SU(2)l quantum numbers such as H^ino and Higgsino. Therefore the 
non-thermal component can explain the dark matter of the universe for m^ ~ 10 3 ~ 5 TeV 
depending on m x . 

For concreteness, let us focus on the case of the H^ino LSP. The annihilation cross 
section is [g3[ (cf. the extensive list |41|| ) 

,2 i3/2 



<? 2 4 1_ 

4 



i 



Mr 



\ a W°W°-*W+W- v i ~ 2vrm 2 r " n2 ' t2l) 



Mr 



In Fig. Ill we show the relic Wino density O x h 2 (where h ~ 0.7 is the present normalized 
Hubble expansion rate |3(|) as a function of m^,. To produce Fig. H, we solve the set 
of Boltzmann equations (|l~8|) (extended to include the charged H^ino NLSP) and take 
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Figure 1: Q x h 2 as a function of the mass m^ for various TFlno masses m x . We use 
for the mass difference between charged and neutral IFinos Am x = 165 MeV, and 
fix M = Mp. The shaded bar corresponds to the 2<r region of Ocdm as reported by 



WMAP [36]. As usual, we fix the present energy density of the universe so that f2 can 
formally become larger than one ('overclosure'). 



into account coannihilation. The coannihilation effect between LSP and NLSP becomes 
important for large m^, and explains the deviation of the contours from straight lines 
(see in particular the m x = 100 GeV contour in Fig. |l|). If Td > m x /30, the VFinos are 
in thermal equilibrium, and the relic abundance does not depend on tha, any more. This 
explains why the contours become horizontal for large m^ in Fig. |l]|j For Td < m x /30, the 
estimate ( p0| ) turns out to be a reasonable approximation. In particular, the temperature 
Td can be as low as 100 MeV without overdosing the universe. If we fix M to be the 
Planck scale, we find that m^ should exceed 10 3 TeV. 

It is instructive to study the dependence of the baryon density in Eq. (|l7]) and the 
relic dark matter density in Eq. ( Jl9| ) on the physical parameters. For case A in Eq. (|17[), 
amazingly, the ratio is independent of the <fi mass, 



»v 



k\ U(av)£m 3/2 



i ( MY 
\~Mp~) 



10 4 \k\ 



m. 



/ M 
Imp 



(22) 



n b \ "i / \MpJ V TO 3/2. 

Our scenario relies on a large gravitino mass which is realized in anomaly mediation 



|1C, 11] where the Wino mass is suppressed by a loop factor, e.g. m x /m 3 / 2 ~ ^/(I^tt 2 ) 
for the WTno \- This implies n x /n\ ) ~ (few) x 10~ 2 for M ~ Mp. Hence, our scenario 



3 For such large m^, the standard WIMP scenario works and thus the Bino is a good dark matter candidate. 
However, this would require even larger <f> masses, and according to Eq. (fiT\), large m^ would make baryogenesis 
difficult. 



predicts for the ratio of dark matter to baryon densities 

-^ ~ |k| _1 x few x 10" 2 x * — x . (23) 

\l b m nuc l C on \M P J 

In particular, for m x of the order 100 GeV (and M ~ Mp), the observed ratio OcDM/^b — 
5 [|36| finds a very natural explanation within the framework described here. The same 
is true for the Higgsino LSP case since it naturally has a mass of the order of the lyino 
mass. 

It is interesting to relax the assumption M ~ Mp and take, for instance, M to be of 
order GUT or compactification scale, M ~ Mgut — 3 x 10 16 GeV, or the string scale. If 
so, the (ft mass can be substantially lower, m^ ~ m 3 / 2 ~ 100 TeV. 

Concerning cold dark matter, our analysis coincides with the one of Ref. |jj|] if we 
identify (ft as a modulus. However, our assumption that (ft is odd under .R-parity completely 
avoids the gravitino problem caused by the (ft decay [42, 43 1. 



4 Conclusions 

We have discussed a scenario where both the observed baryon asymmetry and the cold 
dark matter originate from decays of a heavy scalar field (ft. In an example we imposed Zg 
symmetry such that (ft number is effectively conserved in the (ft oscillation era. This allows 
to have a initial asymmetry q^ to be conserved until the (ft decay, and gu is converted into 
baryon asymmetry. The baryon asymmetry is automatically in the right ballpark for (ft 
masses which are high enough to evade the BBN constraints. The (ft decays also produce 
LSPs. For sufficiently high decay temperature T^, pair annihilation is still partially effec- 
tive if the LSP is the Wino or the Higgsino, and hence the number density of LSPs gets 
reduced into a correct size for the dark matter, n x <C n&. In one of our scenarios (case 
A) we require a heavy gravitino, m 3/ / 2 ~ 100 TeV, loop contributions to gaugino masses 
become important, so that one naturally obtains a lyino or Higgsino LSP. 

Amazingly, in that context, our mechanism predicts Q x /Qf, ~ few x 10 -2 x 
(m x /m nuc i con ) independently of the (ft mass m^. It can hence very naturally account 
for U x ~ 5 Qf, although each (ft decay produces (at least) one superpartner and ttjlsp > 
100 x m nuc i eon . In our scenario, the role of the (ft field is thus twofold: it serves both as 
source of the observed baryon asymmetry and explains why H^inos or Higgsinos can be 
cold dark matter despite of the large annihilation cross section. 

Our mechanism also opens new possibilities in inflation model building. One could, for 
instance, envisage a sequestered scenario where the inflaton is (geometrically) separated 
from the MSSM fields. This is usually a problem because then the inflaton reheats the 
hidden sector. In our scenario, however, such a setup may be viable because (dark) matter 
and radiation are generated by the (ft field rather than inflaton decays. 
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